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ABSTRACT 

Recent observations of partial X-ray eclipses from 4U1822-37 have shown 
that the central X-ray source In this system Is diffused by a large 
Compton-thick accretion disk corona (ADC). We show that another binary, 
4U2129+47, also displays a partial eclipse and contains an ADC. The possible 
origin of an ADC Is discussed and a simple hydrostatic evaporated ADC model Is 
developed which, when applied to 4U1822-37, 4U2129+47 and Cyg X-3, can explain 
their temporal and spectral properties. The quasl-sinusoldal modulation of 
all three sources can be reconciled with the partial occultatlon of the ADC by 
a bulge at the edge of the accretion disk which Is caused by the Inflowing 
material. The height of this bulge Is an order of magnitude larger than the 
hydrostatic disk height and Is the result of turbulence In the outer region of 
the disk. The spectral properties of 4U1022-37 and Cyg X-3 are very similar, 
while 4U2129+47 exhibits a dissimilar spectrum which appears characterl stIc of 
Sco X-l/galact1c bulge sources; the spectra of all three sources can be 
reconciled with Compton scattering of the original source spectrum by the 
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AOC. Spectral variations with the overall Intensity state of Cyg 
probably caused by changes In the optical depth of the corona. A 
of our model Is that any accreting neutron star X-ray source In a 
semi-detached binary system which Is close to Its Eddington limit 
contain an optically thick ADC. 

Subject headings : stars: accretion - X-rays: binaries 
Dept. Physics and Astronomy, Univ. of Maryland 
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I. INTRODUCTION 

The 5.57-hr X-ray binary system 4U1822-37 (Mason et a1. 1980) displays a 
very unusual partial X-ray eclipse, with 50% of the X-ray flux remaining at 
minimum light (White et al. 1981; hereafter Paper 1). The energy independence 
of this eclipse led us to conclude in Paper 1 that the compact X-ray source is 
diffused by a large (~ 0.3 R^) Compton-thick highly ionized corona, above and 
below the center of its accretion disk. This was not the first suggestion for 
accretion disk coronae (ADC) in compact X-ray binary systems. Jones and 
Foreman (1976) proposed that the residual emission from Her X-1 during the off 
part of its 35-d cycle is caused by an ADC scattering X-ray emission from the 
pulsar around the occulting twisted disk. Subsequent observations have 
supported this hypothesis (Becker et al. 1977; Pravdo et al. 1979; Parmar, 
Sanford and Fabian 1980; Bai 1980). In Cyg X-1, Liang and Price (1977) and 
Bisnovatyi -Kogan and Blinnikov (1977) reconcile the hard X-ray spectrum of 
this source with disk accretion onto a black hole, by invoking the 
upscattering of a central source of soft photons in a 10^ K ADC. These 
earlier results provided only implicit evidence for the presence of ADCs, 
whereas 4U1822-37 furnishes a much less subjective experimental signature and 
indicates that ADCs may have a significant effect on our current understanding 
of binary X-ray sources in general. 

He will show that another X-ray source, 4U2129+47, is diffused by an ADC 
and hence appears extended to an outside observer (II). A simple model for an 
evaporated ADC is developed and shown to be consistent with the observed 
coronae (III and IX). The observed eclipse duration and depth are used to 
determine the coronal radius and the nature of the secondary in 4U21 29+47 and 
4U1822-37 (IV). 

In addition to the eclipse, 4U1822-37 also exhibits a 50% peak-to-peak 
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secondary modulation, with its minimum occuring at phase ~ 0.2 prior to the 
eclipse center. In paper i we pointed out that this modulation is very 
similar to the 4.8-hr modulation of Cyg X-3. A similar effect is also 
apparent in new X-ray observations of 4U21 29+47. We will further explore the 
possibility, discussed in Paper 1, that the secondary modulation results from 
the occultation of the ADC by structure at the periphery of the accretion disk 
caused by the inflowing gas 'cr-jm (V and VI). The spectral properties of Cyg 
X-3, 4U1822-37, and 4U2129+47 are presented and discussed in the context of 
previous models for Cyg X-3 and the ADC model proposed here (VII). Finally, 
the implications of our results as regards the nature of the galactic bulge 
sources is discussed (X). 

II. NEW 4U21 29+47 DATA 

The optical light curve of 4U2129+47 is modulated with a period of 5.2 
hrs and an optical spectrum which is cooler at minimum light suggests that it 
results from the rotation of an X-ray heated secondary, similar to HZ Her 
(Thorstensen et al. 1979). This is in contrast to the optical counterpart of 
41)1822-37, which shows the color-independent eclipse of an accretion disk 
(Mason et al. 1980). Ulmer et al. (1980) have reported evidence for an X-ray 
eclipse of 4U2129+47 coincident with the optical minimum. Using the HEAO-1 
A:"*^ experiment we have obtained three 10-day scanning observations in Dec 


■•■The A2 experiment on HEAO-1 is a collaborative effort led by E. Boldt of 6SFC 
and G. Garmire of CIT, with collaborators at GSFC, CIT, JPL and UCB. 


1977, June 1978 and Dec 1978 and one 3-hr pointed observation on 1978 Nov 28 
(see paper 1 for instrument details). In May 1979 we made a 4-hr observation 
using the SSS and MPC detectors on the Einstein Observatory HEAO-2. 
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The folded HEAO-1 A2 count-rate data are shown In Figure 1, and the 
Einstein data In Figure 2 , both with the non-source background removed, and 
plotted against the orbital period ephemerls given by Thorstensen et a1. The 
X-ray eclipses at phase zero are evident In all the data, with a 251 residual 
flux at minimum. The eclipse data with the best statistical quality that we 
have, from the MPC/SSS, show an extended Ingress lasting at least 0.1 In phase 
(Figure 2). The MPC hardness ratio (the flux In the 3.3-10.0 keV band to that 
over 1.1-3. 3 keV) Indicates a slight hardening of the spectrum through the 
eclipse. Spectral fits to the MPC data show a small Increase In absorption 
from <1.3 X 10^^ to 3.2 i 1.0 x 10^^ H atoms cm"^, while similar fits to the 
SSS data give only an upper limit to of < 2 x 10^^ H cm”^. 

Three possible explanations for a partial eclipse were discussed In paper 
1: a) The grazing occultatlon of a stellar X-ray source by Its companion, b) 

Immersion of the system In a stellar wind, or c) an ADC surrounding and 
diffusing the X-ray emission from the source. The energy Independence of the 
partial eclipse of 4U1822-37 eliminated the first two models for that 
system. Even If the small Increase In absorption seen by the MPC from 
4U2129+47 during eclipse Is significant. It Is still at least one order of 
magnitude below that expected from the first two models (see paper 1), so we 
conclude that the emission from 4U2129+47 Is also diffused by an ADC. 

III. ACCRETION DISK CORONAE (ADC) 

Several workers have proposed that the accretion disk surrounding a 
compact X-ray source will be sandwiched between a hot corona. Models for 
generating an ADC can he divided Into two groups: the dissipation of acoustic 

or magnetic energy flux generated by convective or viscous turbulence within 
the disk, or the evaporation of disk materia) by the central source. 

Icke (1976) first suggested that tlie viscous turbulence within an 
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accretion disk, responsible for the transport of angular momentum outwards, 
could also give rise to acoustic noise In the disk and generate a hot corona 
(analogous to similar models for the solar case). In the same vein, 
Bisnovatyl-Kogan and Blinnikov (1977) and Liang and Price (1977) discuss 
convective turbulence In the inner photon-pressure-dominated region of an 
accretion disk (see also Galeev, Rosner and Valana 1979), although the photon 
pressure region Is In the Innermost region of a disk and Is much smaller than 
the coronal radii under consideration here. Most of the accretion disk 
surrounding any compact object (Including a neutron star or white dwarf) will 
be gas-pressure-dominated, and It Is unclear whether convective Instabilities 
can generate a corona In such a situation (Tayler 1980; Shakura, Sunyaev and 
Z11 Itinkevich 1978; Livio and Shaviv 1977), or what the contribution from 
viscous turbulence will be (Icke 1976). In addition, there may be magnetic 
loops threading the disk, again analogous to the solar ca^e. Because of these 
uncertainties we will concentrate on the other possibility, an evaporated ADC. 

Shakura and Sunyaev (1973) pointed out that the central source in an 
X-ray binary system will evaporate material from the surface of its accretion 
disk. If It does not escape from the system, this hot material then settles 
above and below the disk as a corona. The physics of this process has been 
thoroughly discussed with regard to the self-excited stellar wind mechanism 
for causino mass transfer in X-ray binaries (e.g. London, McCray and Auer 1981 
and refs therein). These calculations can be used to predict the properties 
of an evaporated ADC. 

When the atmosphere of a star or an accretion disk is illuminated by a 
source of X-rays it will photo-absorb the incoming radiation and re-emit the 
energy as continuum and line emission. The brenisslrahlung cooling per unit 
mass is proportional to n/T, while the heating per unit mass depends only on 
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the Ionization state of the gas, not Its density; therefore n « T ' . The 
equation of state, however, gives n « PT"^, I.e. n « P"^ and, for hydrostatic 
equilibrium, the gas Is thermally unstable. At sufficiently low densities an 
equilibrium can only be achieved by abruptly raising the gas temperature to 
change Its Ionic state, thereby reducing the X-ray heating (e.g. McCray and 
Hatchett 1975). The presence of line cooling will complicate this simple 
picture somewhat, but does not change the overall result. If the thermal 
energy of the gas exceeds the gravitational energy then a stellar wind Is 
formed. If not, surrounding an accretion disk there may be a large volume 
where material within the radius - IMj^/Mq). T^^ Rq will remain In 
hydrostatic equilibrium. The coronal temperature Ty In units of 10^ K Is 
given by the balance between bremsstrahlung losses, Compton cooling and 
Compton heating. For the case in which we are Interested, the bremsstrahlung 
Is negligible and, from Shakura and Sunyaev (1973), we obtain Ty = 0.25 T^, 
where T^ Is the central source temperature. T^ Is unknown but we will as:ume 
It to be similar to that measured from X-ray pulsars and set T^ = 24 
(appropriate for Her X-1), consistent with the value used by McCray and 
Hatchett, but note that the value of T^ may noticeably change the size of the 
coronal region from source to source. The density structure of the corona as 
a function of the angle from the plane of the disk 0 at a radius R can be 
obtained by equating the gas pressure to the net gravitational and radiation 
pressure (after Lynden-Bell 1969). The resultant gas density n is 
Gausslan-dlstributed about the accretion disk plane as, 



where 
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and Lg^jj Is the Eddington-limited luminosity for a source of mass Mj^. The 
critical density at which line cooling becomes Important Is the maximum 
density n^^^ possible at the base of the corona where the high temperjuture Is 
maintained, n^^^ Is determined solely by and R, the luminosity of and the 
distance from, the X-ray source. This corresponds to an Ionization 
parameter c of ~ 10^ (McCray and Hatchett 1975) such that 
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kT^R“ 


for a 1 .'I Mq neutron star. 

An estimate of the optical depth t through the corona at high 
Inclinations (1=90°- 0) may be obtained by Integrating (2). The luminosity 
Lj., above which the corona becomes optically thick Is 

l-c • '■edd • '"(n'",' "-9 s-‘ (3) 

where R^^ Is the Inner radius of the corona. The low density of the corona 
will probably prevent It extending all the way In to the Alfven surface and It 
will be disrupted by the magnetic field at some point within the "transition 
zone" (discussed In Ghosh and Lamb 1980) and R^p will be ~ 10"^ R^. Therefore 
an ADC will be Compton-thick only for i 0.1. The optical depth as a 

function of Inclination Is shown In Figure 3 for an arbitrary value of R^p * 
0.05 Rq and Ty of 6. For L^/Lg^j^ > 0.1 the source will always be 
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obscured (t > 1) above Inclinations of 60°, and might possibly be at lower 
Inclinations as well, depending on the density (which Is specified by L^). 

For T > 1 the photon density at a given distance R will be enhanced because 
photons will be trapped In the corona for the time that It takes them to 
escape. Ross (1979) has shown that In this case the true luminosity of the 
X-ray source cannot be ascertained without concurrent knowledge of the density 
structure. For our discussion the Interior modelling of the kdC Is not 
important; we require only the conditions to produce an ADC. In Figure 4 the 
probable appearance of the evaporated ADC corona of 4U1822-37 Is Illustrated. 

IV. THE ECLIPSE BY THE SECONDARY 

All the eclipse :alculat1ons have been performed via numerical rather 
than analytic means. In order to model the eclipses observed In these systems 
It Is useful to define the dimensions of the optically thick region of the 
corona as they would appear to an outside observer. This can be done by 
considering the "surface of last scattering" defined as the radius R^ 
(perpendicular to the line of sight) and height H^. (parallel to the disk) at 
which T » 1 (cf Fabian and Pringle 197B). Numerical Integration of (1) and 
(2) gives 



(4) 


Most of the subsequent analysis Is relatively Insensitive to uncertainties in 
(4) but we shall Indicate where our results depend on It. Two further 
assumptions were made: (1) a spherical secondary that fills Its Roche lobe, 

so that Its mass and radius Rj are related by the Roche mass-radius 
relations giver by Paczynski (1971), and (2) the x-ray source Is an accreting 
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neutron star with 0.8 < < 2.0 M^. These two assumptions, Keplers third law 

and (4) make the duration and depth of the eclipse functions of only 1, and 

Re* 

The results of fittinq the secondary eclipses of 4U1822-37 (taken from 
the HEAO-1 A2 pointed observation, paper 1) and 4U2129M7 (MPC data, Flqure 2) 
are given in Figure 5 as Rj vs 1 and R^ vs 1. From Figure 5 It Is apparent 
that the radius of the ADC corona Is Insensitive to both the companion radius 
and the system Inclination, and that the coronal radius in 4U1822-37 at least 
half that of 4U2129+47. The constraints imposed by the eclipse on the nature 
of the secondary 1^ these two systems are as follows: 

(1) 4U1822-37: The statistical errors are dominated by the uncertainty 

In the mass of the compact object (assumed in Figure 5 to be 0.8 < < 2.0 

Mq). The error associated with any possible uncertainty In the coronal height 
Is much smaller than this. A hydrogen zero age main sequence (ZAMS) secondary 
gives 1 - 75°. If the Inclination Is > 85° then the mass of the secondary 
will be less than 0.085 M^, tne minimum value for hydrogen nuclear burning. 
However, a degenerate secondary that fills Its Roche radius exceeds the 
maximum theoretical radius for a cold hydrogen white dwarf ( ~ 0.1 Rq, 

Zapolsky and Salpeter 1969) and limits 1 < 85°. For 1 < 67° the secondary 
would be a helium star, however the detection of hydrogen emission In the 
optical spectrum makes this unlikely (Mason, private communication). 

(2) 4U2129+47: We have high quality (pointed) coverage for only one 

third of the eclipse and the statistical errors for R^. and Rj are comparable 
to those from the uncertainty In the mass of the X-ray source. The error- 
region In Figure 5 for this source Is defined from only the statistical 
uncertainties, and the error region moves to the left or right as the X-ray 
source mass changes. Because the coronal radius Is so large this solution Is 
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sensitive to the assumed coronal height. We used a height of 0.4 R^; for 
lower values the solution does not change significantly, but for greater 
values It gives larger secondary radii for a given Inclination, degenerate 
secondaries are ruled out for all possible solutions. A heilum star Is also 
unlikely because of the observed color variations and H line emission 
(Thorstensen et al. 1979). A hydrogen ZAMS sequence secondary gives 
1 ~ 82®. A larger coronal height than that assumed will, for a g<ven mass 
radius relation for the secondary, move 1 closer to 90°. We note that the 
shape of the ecripse center Is quite dependent on the coronal height used, so 
that a complete observation of the whole eclipse with high statistical 
precision sh.'jld remove any uncertainty concerning the coronal height. 

V. THE ECLIPSE BY AN ACCRETION DISK BULGE 
The periodic hump seen In the light curves of many cataclysmic variables 
is :)elieved to come from a hot spot at the point of impact of the inflowing 
gas stream with the edge of the accretion disk (e.g. Robinson 1976). Intense 
flickering of the light from the hot spot suggests turbulence ind this may 
cause a thickening of the disk in the vicinity of the hot spot (e.g. Lubow and 
Shu 1976). Chester (1979) and Alpar (1979) have modelled the time delay of 
the optical pulse from OQ Her (outside eclipse) in terms of a reprocessed 
pulse on the outer disk surface and both found a thickening of the trailing 
edge of the disk, although the vertical extent of this structure was not 
determinable. Hassall et al. (1981) have suggested a similar scheme to 
account for the reprocessed optical pulses from H2252-035. Periodic 
absorption dips from 4U1915-05 have been explained In terms of the occultatlon 
of the X-ray source by a structure in the outer edge of the disk (White and 
Swank 1982; Walter et al . 1982). A gradual decline In the light from the 
optical counterpart to 4U1822-37, prior to the eclipse of the accretion disk 
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by the secondary, may be caused by a bulge, near the hot spot, occulting part 
of the disk (Mason et al. 1980), where the hot spot In this case Is relatively 
dark compared to the X-’^ay-lllumlnated disk. This bulge may also occult the 
ADC and cause the sinusoldal-llke modulation of the X-ray flux (Paper 1). 

A number of similarities between the X-ray properties of 4U182Z-37 and 
Cyg X-3 were given In Paper 1, where It was suggested that they might te 
Intrinsically similar systems (see also VIII). Our simple description "a 
quisl-slnusoldal modul^tloi." glosses over several Interesting details that 
deserre more attention, and a seif-consi stent picture should be capable of 
explaining their differences, as well as their similarities. The minimum In 
the light curve of Cyg X-3 Is asymmetric, with a steeper fall than recovery 
(e.g. Mason et al. 1976a, and Figure 9). 4U1822-37 displays a different 

behavior with a more gradual decline to minimum, a steeper rise to a plateau 
where the eclipse by the secondary occurs and then a second Increase to 
maximum light (Figure 6 this paper; Figure 1 paper 1). These features have 
remained constant over many years (Eisner et al . 1980 and paper 1). 

In this section we Investigate the possibility that both these 

modulations are caused by occultat(on of the ADC by azimuthal variations In 

the height of the outer edge of the accretion disk. The numerical code used 

to model the eclipse by the secondary was modified to Include an occultation 

by material at the periphery of the disk. In all these calculations we 

1 125 

assumed the sianoard thin disk radius-height relation of » 0.05 R^* 
(corresponding to an outer disk temperature Tq » 2 x 10^ K) and that the bulge 
Is confined to the edge of the disk. We used the thin disk relation for 
convenience and note that the disk may be thicker. The following results are 
Independent of reasonable variations In the height-radius relation of the 
disk, or of the coronal radius-height relation. 
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A simple bulge profile (shown at the bottom of Figure 6) that had a sharp 
rise at the confluence of the gas stream with the disk and decayed linearly 
downstream of that point was tried first. Two parameters modified the shape 
of the resulting light curve for this geometry: a) the rate of decay of the 
bulge with azimuthal angle and b) the ratio of the radius of the optically 
thick region of the corona to that of the disk. Figure 6 shows the effect of 
Increasing R^., keeping the other parameters, In particular the disk radius, 
constant. The secondary eclipse has been removed for clarity and the HEAO-1 
A2 pointed observation of 4U1822-37 Is given (but not fit) for comparison. 
well Inside R^j gives a sharp minimum, but as the coronal radius apprcsches the 
disk radius the minimum beccr.«o ehallow, with a slight asymmetry similar to 
that of the Lyq X-3 light curve. This simple geometry cannot, however, 
reproduce either the plateau around the secondary eclipse In 4U1822-37 or the 
duty cycle of the modulation of Cyg X-3. 

The plateau In the light curve of 4U1822-37 suggests a thickening of the 
disk near ® = 0, i.e. upstream from tht confluence. Including such a 
thickening in the prebulge region gives light curves such as those shown In 
Figure 7. As R^. approaches R^j this plateau becomes smeared out and again we 
obtain a smooth light curve, but now with a longer duty cycle similar to that 
of Cyg X-3. Figure 7 Indicates that all the gross features of the light 
curves of both Cyg X-3 and 4U1822-37 can be reproduced with the same model, 
with their differences arising primarily from the distance between the bulge 
and the optically thick edge of the corona . 
a) The X-ray light cui^ve of 4U 1822-37 

Formal fits to the data were made with the the angle of confluence 6, the 
bulge height above the plane of the orbit Hg, the phase angle through which It 
linear’y decays 0p, the phase angle of the start of the prebulge 6p and the 


14 


prebulge height Hp as free parameters. The top panel of Figure 8 shows an 

optimized light curve for this geometry. We note that the minimum thickness 

at phase ~ 0.25 Is similar to that found by Chester and Alpar for DQ Her 

(at t ~ 0*7 using their phase convention). However, there are still 

systematic deviations from the predicted light curve that suggest the prebulge 

region Is not constant In height, but decays In the same sense as does the 

primary bulge (with a phase angle through which It decays of 9 q 2)* This 

modification Improves the predicted light curve considerably, as shown In the 
2 

second panel, x formally unacceptable at 410 for 163 dof for this fit, 
but because of turbulence the bulge shape Is probably fluctuating on 
timescales much shorter than an orbital period. A sharp rise was used for the 
start of the bulge at the confluence. In reality this feature must be less 
abrupt and the Inclusion of a more gradual rise did not significantly change 
the resulting light curves. 

It might be argued that the light curve is being over-modelled, i.e. 
given enough parameters any random light curve can be reproduced, but not all 
physically reasonable alternatives lead to better fits. A large number of 
bulge profiles were tried, including gaussian and exponential decays, 
prebulges that Increase In height downstream, and a highly elliptical disk 
(widest perpendicular to the line joining the two stars) but none synthesized 
the features of the observed light curve. 

The only other high quality data of 4U1822-37 available comes from an 
Einstein ff>C observation made one year later. The data are shown in the third 
trace of Figure 8, with the best fit HEAO-1 A2 light curve. The data are 
sparse, but the general features of the light curve are confirmed, although 
there Is evidence that the angl<; 6 may have changed slightly. The best fit 
parameters for 4U1822-37 are given In Table 1. The Inclination was varied 
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from 7f)° to 80°, which qave a range of bulge height of 0.15-0.3 Rq, but no 
other parameters were affected by changing the Inclination. Because the 
coronal radius Is fixed by the secondary eclipse, we can determine a radius 
for the disk of 0.6-0. 7 R^. The proposed geometry for 4U1 822-37 Is 
Illustrated In Figure 4. 

b) The Optical and UV Light Curves of 4U1822-37 

The optical light curve consists of a gradual decline from ~ 0.6 - 
0.9, followed by a relatively sharp eclipse of the accretion disk by the 
secondary. Mason et a1 . (1980) suggested that the pre-eclipse decline Is 
caused by an occultatlon of the central region of the disk by a bulge at the 
edge of the disk. The numerical code was modified to model the expected 
optical eclipse using the azimuthal bulge profile from the X-ray data, 
assuming that the disk Interior (Including the bulge) Is uniformly bright, and 
that the outside rim of the disk, not exposed to the X-ray source. Is dark 
relative to the disk Interior. The optical data from Mason et al . (1980) are 
shown in Figure 8. The dashed line Is the predicted light curve assuming no 
additional light In the system. There is a minimum at phase 0.8, but It Is 
much deeper than the observed decline. This Is because the light curve is 
made up ot a combination of both the occultatlon of the disk by the bulge, and 
the self occultatlon of the bulge itself which, at such high Inclinations, 
contributes a sizable fraction of the light. The effect of adding additional 
light from the secondary Is shown as the solid line; 40% of the total light Is 
assumed to be from an X-ray heated secondary and Is modulated by 100%. The 
resultant light curve Is closer to the predicted light curve, but It still Is 
not satisfactory. Therefore, while the pre-eclipse decline can probably be 
accounted for by the effects of viewing the bulge at different binary phases, 
the white light data must be complicated by either a non-uniform distribution 
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of light across the disk, or additional light from the hot spot/gas stream. 

We can obtain limits on the radius of the accretion disk from the white light 
eclipse In a manner similar to those obtained for the ADC radius, and find 
= 0. 7+0.1, consistent with the value found from the X-ray light curve. It Is 
Interesting to note that a second minimum at phase 0.5 may correspond to a 
partial eclipse of the secondary by the accretion disk (cf. Longmore et al. 
1981). 

Mason and Cordova (1982) have recently obtained a coarsely binned 
measurement of the far UV light curve (in seven binary phase Intervals) and 
find that It resembles the X-ray light curve with a minimum at ♦ ~ 0.9, 
although with a larger amplitude. This Is In good agreement with the 
modelling In Figure 8 and suggests that the far UV comes entirely from the 
X-ray illuminated accretion disk and Is uncontaminated by the other light 
sources. 

In paper 1 evidence was given for a phase shift between the X-ray and 
optical minima, with the optical leading the X-ray by 0.04 In phase. In order 
to correctly center the predicted secondary eclipse, the optical data shown In 
Figure 8 Is phase shifted by 0.015 from that of Mason et al (1980), which is 
within the 0.018 phase uncertainty of the current ephemerls (Paper 1). The 
remaining phase difference between the X-ray and optical epochs indicates that 
the distortion of the eclipse caused by the pre-eclipse decline has biased the 
definition of the center of the optical minimum. 
c) Cyg X -3 

The azimuthal profile for the disk used to fit 4U1822-37 was applied to 
an average light curve of Cyg X-3, taken from MPC data obtained during an 
Einstein observation In Dec 1978 (the same data used by Eisner et al. 1980). 
The source was In a very low state (~ 50 UFU) and exhibited variations In the 


17 


light curve from cycle to cycle (see Figure 1 of Eisner et al.). The data and 
best fit light curve are shown In Figure 9. The general trends of the 

2 

modulation are reproduced remarkably well and the formally unacceptable x 
understandable In terms of variations In the light curve on short timescales 
well In excess of the excellent statistics. Apart from the bulge height, the 
only changes made to the 41)1822-37 azimuthal profile were a decrease In 
0JJ to 100° and an Increase In 6p to 290®. The ratio of the prebulge height 
to the primary bulge was the same for both source fits. Because we have no 
prior knowledge of the coronal radius, we can obtain only a range of bulge 
heights as a function of the ratio R^/Rd which, as expected, was larger than 
the values for 4U1822-37 (~ 0.75 vs 0.5). Figure 10 displays the predicted 
bulge height as a function of Inclination for different disk radii. The 
arbitrary phase In Figure 9 has 4=0 corresponding to approximately where the 
eclipse by the secondary would occur If 6 ~ 40®. 

The lack of any eel lose places an upper limit on the secondary radius as 
a function of Inclination. This Is Indicated In Figure 10 assuming = 1.4 
Mq. The 4.8-hr orbital period of Cyg X-3 means that a Roche-flll Ing cold 100% 
hydrogen white dwarf Is just possible (see Eisner et al . 1980), but the white 
dwarf would be of an extremely low mass ~ 0.005 Mq and on the verge of 
becoming a planet (Zapolsky and Sal peter 1969). Thf Improbability of such a 
large mass transfer rate (L^ ~ 10^® ergs s"M occuring In such a system 
suggests the secondary Is not degenerate, which limits 1 < 80®. A hydrogen 
ZAMS secondary gives 1 < 68®. 

The average luminosity of Cyg X-3 in the 2-10 keV band can vary by a 
factor 5 on timescales of months. Five days of HEAO-1 A2 scanning 
observations of Cyg X-3 In Its high state In 1978, folded at 4.8-hrs, are 
displayed In Figure 9. The scatter In the data points Is greater than that In 
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the MPC data because the observation is over such a long time interval and the 
source is only seen for ~ 1 minute out of every 40, so that it is very 
susceptable to local fluctuations in the source intensity. Nonetheless the 
general shape demonstrates that there have been no major change in the coronal 
radius or the bulge height. 

Determination of the bulge height above the disk in Cyg X-3 is dependent 
on knowledge of the disk radius and system inclination. From Figure 10, the 
peak bulge height in Cyg X-3 must be at least half of the disk radius or > 15 
times the hydrostatic disk thickness (for Tg = (.6-2) x 10^ K) and subtends an 
angle of > 40° to the central source. This compares with the 4U1822-37 value 
of about one third the disk radius or ~ 7-15 times the hydrostatic thickness, 
subtending an angle of ~ 20° to the central source. 
d) 4IJ2129-^47 

The statistical quality of these data do not constrain the model as 
severely as do those from Cyg X-3 and 4U1822-37, so that no rigorous attempt 
has been made to model its light curve. However, the overall light curve of 
4U2129+47 seen on 1977 day 342-352 does bear a marked resemblance to the 
quasi-sinusoidal modulation of 4U1822-37, with a second minimum occuring at 
phase ~ 0.65 (Figure 1). This minimum is not as obvious six months later 
(where the statistics are worse), and is not at all evident in either the 
pointed or scanning HEAO-1 A2 data when the source was at a factor 2 lower 
intensity. In May 1979, the source returned to its higher intensity and a 
second minimum again seems to be present, but at a slightly later phase 
of ~ 0.75 (Figure 2). The possible shift in the phase of the minimum 
indicates the tantalizing possibility of major changes in the outer disk 
structure, possibly related to the change in X-ray intensity. This source 
clearly deserves further attention from X-ray observers. 
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VI. THE NATURE OF THE BULGE 

The accretion disks In systems such as these are a consequence of the 
necessity for the loss of angular momentum before accretion onto the compact 
object can occur. Flannery (1976) and Lubow and Shu (1975) suggest that the 
equilibrium outer radius of the disk R|j| Is given by the point at which the 
specific angular momentum of the disk material equals that of the Infalling 
material. If the disk radius differs from this, then a large velocity 
differential occurs across the confluence and causes strong turbulence that 
tends to spin the disk up or down, accordingly. Tabulated values for both Rj^ 
and the angle of confluence 6 are given In Lubow and Shu and Flannery and we 
can compare them with the observed disk values for 4U1822-37. Mass ratios 
in the range 1 to 10 give Rm between 0.15 and 0.49 Rq, which is 
significantly less than the disk radius estimates for 4U1822-37 of 0.6-0. 7 
Rq. Further evidence that the disk in 4U1822-37 Is larger than the 
equilibrium value is given by 6 which is ~ 20° below that expected for the 
given mass ratio range. So, unless the mass ratio in this system is extremely 
large (i.e. the compact object is '' 2 M^), the disk radius must be about 
twice R|v|. Therefore, large slip velocities and turbulence may be expected at 
the confluence. 

Turbulence will increase the effective gas pressure within the disk and 

cause the disk to thicken at its edge. The turbulent pulsation gas pressure 
2 

is given by pV^ where p is the density and the mean velocity of the 
turbulent motion. This Increases the scale height of the disk by a factor T^ 
given by 



■ U) 


for u) > 1 


( 6 ) 
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where tu 1s the turbulent Mach number and y the ratio of specific heats 
(Shakura, Sunyaev and Zllltinkevich 1978). Is unknown but should be 
comparable to the velocity difference aV across the confluence and must be 
less than the keplerlan velocity. The sound speed In the disk will be ~ 10 km 
s"^ and Is much less than the outer disk keplerlan velocity of ~ 600 km s"^. 

An (i» of 7 to 15 for 4U1822-37 gives a ~ 100 km s"^. The factor-of-two 

larger bulge height for Cyg X-3 can be accounted for by a different mass ratio 
and/or disk size, either of which can Increase aV across the confluence. At 
such large Mach numbers the gas pressure will be Increased by vorticity caused 
by the generation of shock waves In the turbulance, further thickening the 
disk. It Is unclear how far Into the disk the turbulence will penetrate. 

Such large supersonic turbulent velocities must lead to fluctuations In the 
size and shape of the bulge consistent with the observed cycle-to-cycle 
variations In the structure of the light curve of Cyg X-3. These variations 
tend to occur predominantly on the rise and peak of the 4.8-hour cycle 
(van der Klis, priv. comm.; Parsignault et al . 1976), suggesting that the most 
dramatic variations In the disk height occur In the region upstream from the 
confluence, I.e. In the structure of the second bulge. The longer term 
variations In the light curve seen by van der Klis and Bonnet-Bldaud (1981) 
may reflect changes In the coronal radius, disk radius or disk structure. 
Bonnet-Bldaud and van der Klis (1981) have also recently found evidence for a 
20-day 1% periodic phuse shift In the times of minimum light from Cyg X-3. 

This suggests that there may be a slight oscillation In the position of the 
bul ge. 

The accreting material Is spun up to the disk velocity and Injected Into 
orbit about the compact star, but the second bulge Indicates that It Is an 
orbit that Intersects stable orbits of disk material, which causes It to 
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collide with the disk a second time, giving rise to another shock (similar to 
shock II In Lubow and Shu). The three-body hydro-dynamics of this situation 
are complicated, however, we note that this first direct measure of the 
azimuthal structure at the edge of an accretion disk provides an Important 
constraint to any general description of disk accretion. 

Smak (1979) and Sulkanen, Brasure and Patterson (1981) have measured the 
radii of the accretion disks In several cataclysmic variables (CVs) and found 
them all to be substantially larger than R|^. We might, therefore, expect 
large bulges In CVs as well, and It Is perhaps a little puzzling why such 
prominent features may have gone untoticed In such well -studied systems. The 
outer thickness of the disk Is determined by Its temperature and part of the 
explanation may come from X-ray heating enhancing the thickness of the outer 
disk structure. Further, for high Inclinations the luminosity of the hot spot 
can be comparable to that of the disk (Bath et al . 1974), and the 
characteristic hump In the light curve from the spot could be masking the 
pre-eclipse decline (like that seen In 4U182Z-37) caused by the 
sel f-occul tatlon of the disk by the bulge. 

RW Tri Is a CV whose orbital period (5.57-hrs) and dynamic properties are 
very similar to those of 4U1822-37. Longmore et al. (1981) have recently 
observed an IR eclipse of the secondary by the accretion disk and from this 
they estimate a radius for the disk of ~ 0.75 Rq, or five times Rf^j. They also 
measure a "surprisingly large" thickness for the disk of 0.11-0.17 Rg. The 
disk thickness we measure In 4U1822-37 at ♦ ~ 0.0 and ♦ ~ 0.5 (see the bulge 
profile In Figure 8) Is very similar to that seen by Longmore et al . 

Therefore, such large bulges may be a common phenomenon In CVs and low mass 
X-ray binaries. 


VII. SPECTRAL RESULTS 
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N'>n-d1spers1ve X-ray spectra of Cyg X-3, 4U1822-37 and 4U2129+47 were 
obtained using HEAO-1 A2 and the Einstein SSS. The HEAO-1 A2 data come from 
two Independent detectors sensitive In overlapping and complementary energy 
bands 2-20 keV and 2-60 keV, each with 64 channels of PHA. The SSS covers the 
range 0.5-4. 5 keV with 128 channels of PHA, and has a factor 4 better energy 
resolving power at 3 keV than do the HEAO-1 A2 detectors. Additional spectral 
Information out to - 20 keV Is provided by eight channels of PHA from the 
Einstein MPC. AH uncertainties are 90% confidence, 
ll Cyg X-3 

The spectrum of Cyg X-3 changes radically from Its low to high state 
(Serlemitsos et a1. 1976; Sanford, Mason and Ives 1976) and this has caused 
some confusion regarding the spect'^al properties of this source (cf Becker et 
a1. 1978; Bllsset, Mason and Culhane 1981) 

The low state HEAO-1 A2 results come from a pointed observation on 1978 
Dec 2. The phase averaged spectrum Is a power law, with an energy Index a of 
0.1 i 0.1 similar to that found by Becker et a1 . (1978). The extended energy 
range of HEAO-1 A2 reveals a break at ~ 18 keV above which energy the spectrum 
steepens with an e- folding energy of - 16 keV. The MED fit requires an 
additional soft excess with a thermal bremsstrahlung temperature of ~ 0.3 keV, 
3 times lower In temperature than that reported by Becker et a1. 11978). This 
excess may come. In part, from stars In the nearby Cyg 0B2 association Icf. 
Hamden ot a1 . 1980). The photo-electric absorption, kept the same for both 
components in the spectral fitting, Is - 7 x 10^^ H cm"^. 

The Iron line EW Is 1340+60 eV (0.014+0.003 ph cm'^ s"M at 6.76t0.03 keV 

with a FWHM of 1. 0+0.1 keV, again in reasonable agreement with that seen 

2 

previously (Becker et a1. 1978). The x for the fits from both detectors are 
poor, with 335 for 36 dof and 931 for 46 dof In the MED and HED, 
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respectively. The residuals indicate this to be the result of poor modelling 

of the iron line and the high energy break. As the iron line used is a simple 

gaussian, more conplicated line blends can be approximated with a sum of two 

2 

lines, one narrow and one broad, which gave x of 10^ aod 330 in the HED and 
HED, respectively. The EWs of the narrow and broad component are 830 ± 40 and 
840 1 50 eV (0.0081 ph cm"^ s"M at energies of ~ 6.78 and 6.51 keV. The 
width of the broad component is 2.4 keV. The deconvolved spectrum, using the 
single line fit, is shown in Figure 11. 

A sensitive model independent indicator of spectral variability across 
the 4.8-hr modulation is the ratio of the rate in each PHA channel obtained at 
maximum, to that obtained at minimum light. These are shown in the left hand 
side of Figure 12 for the HEAO-1 A2 data such that an Increase in the depth of 
modulation causes the ratio to Increase (i.e. more absorption at minimum light 
causes a larger PHA ratio at low energies). The following conclusions may be 
drawn from Figure 12. 

1. The ratio decreases at low energies, suggesting that the soft component 
is modulated differently than is the higher energy component. 

2. There is a dip around -6.7 keV, indicating a difference in the 
modulation of the iron line relative to the local continuum. 

3. The spectrum hardens at minimum light. 

These points are essentially in good agreement with the conclusions of 
Becker et a1. (1978). Because of the Interdependence of the parameters it is 
difficult to quantify these variations, in particular the variations in the 
iron line. For a single gaussian line there is no formal statistical evidence 
for a change across the modulation in either the line energy or the EW. The 
twu-fcmponent line again ;ives the same EW (within " 10%) but there is 
evidence that the centroid energy of the broad component had increased 
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to ~ 6.7 keV at minimum light, similar to the shift seen by Becker et a1. 
M978). 

The 7 arc mm field of view of the SSS (compared to ~ 3° lor HEAO-1 A2), 

did not Include any known sources In the Cyg 0B2 association. The observation 

was, In part, simultaneous with the HEAO-1 A2 measurement just discussed. The 

ratios of the PHA rates at minimum and maximum light, shown In Figure 13, 

again decrease at low energies and suggest the presence of an unmodulated soft 

t^.^ponent. Formal fits to the total data set for the model of a single power 

2 

law plus absorption model give an unacceptable x of 285 for 63 dof, with a 

large fraction of this x associated with the structure In the spectrum near 

~ 2.5 keV (Figure 14). The Inclusion of a soft thermal component with 

2 

kT ~ 0.3 keV gives a x of 238, and the further Inclusion of a broad sulfur 

2 

line at ~ 2.5 keV reduces x To 101. This fit is fhown in the middle trace of 
Figure 13, where the unmodulated soft component replicates the softening In 
the PHA ratio. The EW of the sulfur line relative to the continuum Is 177 t 
30 eV at an energy of 2.54 t 0.04 keV, with a width of 0.2/ t 0.05 keV. The 
luminosity of the soft component Is a factor two below that derived from the 
MED data, and almost an order of magnitude below that required by Becker et 
al. 

Even with such a small FOV there still exists the possibility that the 
soft component Is not Intrinsic to Cyg X-3, but rather comes from a 
contaminating star within 4 arc min (e.g. Figure 1, Hamden et al.). The 
temperature of - 0.3 keV is typical of that seen from OB stars (Swank, priv. 
comm.) and with ~ 50% lower column density than Cyg X-3, the luminosity Is 
similar to the X-ray bright stars seen by Harden et al. The measured 
absorption for Cyg X-3 corresponds to (7t3) x 10^^ H cm”^, where the 
uncertainty Includes the possibility that the soft emission Is not Intrinsic 
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to Cyg X-3. 

High state spectral data were obtained fron, HEAO-1 A2 scanning 
observations made between 1977 Nov 16 and 26. The PHA ratios of the maximum 
to the minimum phases for the high state are also shown In Figure 12, and the 
most obvious difference from the low state Is that the unmodulated s'^ft excess 
Is not apparent In the high state data. The other aspects of the variability 
of the ratio with energy are similar In both states. Including the dip at the 
1rcn line and the tendency for a harder spectrum at minimum. Acceptable fits 
to the data can be obtained with a relatively flat power law plus a steeper 
second component. For reasons discussed In the next section we assume that 
the steep component Is a blackbody. The best fit gives a blackbody 
temperature of 1.2 ± 0.1 keV In both detectors, and a of 0.6 t 0.1 In the MED 
and 0.9 i 0.1 in the HED. The Iron line EW Is 57C t 190 eV (0.039 t 0.011 
cm”2 s"h centered at 6.56 t 0.15 keV, with a width of 1.3 i 0.4 kcV. 

The for 36 dof, and 61 for 17 dof, for the MED and HED, 

respectively. These results agn>rf with those found earlier for the high state 
by Bllsset, Mason and f'ulhane (1981), except that we see no evidence for any 
phase related variations in the blackbody component. The deconvolved spectra 
for both the high and low states are shown In Figure 11. 

SSS observations of the high state were made on 1979 May 9 and 10. A 
blackbody model gives the b<*st simple fit, with kT - 1.5 keV and a of 507 
for 55 dof. The Inclusion ur a second continuum component does not 
substantially Improve the fit. Sulfur and silicon lines at 2.61 and 1.6 keV, 
each with EHs of 40 i 20 eV and line widths of - 0.2keV, do Improve the fit, 
although not so nijch as to make It formally acceptable. The best fit high 
state column density Is (3t1) x 10^^ H cm"^. In reasonable agreement with the 
radio measurement (Chi: and Bleging 1973). The PHA ratio shows no evidence for 
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any Increase (or decrease) at Tow energies. The upper limit to any phase 
related change in absorption is < 1. 5 x 10^^ H cw*^ . This takes into account 
any uncertainty in the underlying continuum. 

For a distance of 10 kpc (Chu and Sieging 1973) the 1-60 keV luminosity 
is 4.4 X 10^^ erg s”^ in the low state and 1.2 x 10^® ergs s"^ in the high 
state . A1 1 of the evidence is consistent with the additional soft component 
observed in the low state being external to the source. Observable in the low 
state only, its apparent lack of 4.8-hour mod'ilation, in all the data 
available to us, suggests that its origin may lie outside Cyg X-3. If it were 
truly associated with Cyg X-3, and it had a blackbody spectrum for which the 
high energy tail is observable in the X-ray band as the soft component, this 
component would dominate the total source output with a luminosity in excess 
of 10^^ erg s"^* blackbody component alone, high-state blackbody 

component is modulated with the 4.8-hour variation and would totally mask the- 
presence of the soft component seen in the low state in the X-ray band. 

These results can be summarized as follows; 

1. Cyg X-3 shows two distinct intensity states. The low state is an a ~ 0 
ptA^er law with high energy cutoff at ~ ’9 keV. In the high state the power 
law is steeper, a ~ I, and there ib a second - 1.2 keV blackbody component. 

2. No variation in absorption across the modulation is observed. 

3. There is an Intense 1-2 keV broad iron line at ~ 6./ keV, with a possible 
slight dependence In the energy centroid of the line feature with the 4.8-hr 
modulation. The EW in the high state is ~ 700 eV and in the low state is 

~ 1600 eV. The number of photons in the line increases from the low to high 
state by a factor 5. 

4. There is a . arrower, ~ 0.3 keV wide, sulfur line feature at - 2.S keV. 
While the number of line photons increases from the low to high state, the EW 
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decreases from 177 t 30 eV to ~ 40 ± 20 eV. 

5. In the low state there is a bright (L^ ~ 10^^ ergs s"^ if truly 
associated with Cyg X-3) unmodulated soft excess (kT ~ 0.25 keV) that is not 
related to the blackbody component in the high state spectrum. 

^ 4U1822-37 

In paper 1 we found that: 

1. This source is stable over many years and shows a flat, a ~ 0, power law 
with a high energy cutoff at - 18 keV. 

2. The secondry modulation is absorption independent and the column density 
is - 3 X 10^* H cm"^. 

3. An 1100 eV EW iron line is present with a FWHM of 4.4 ± 2.0 keY. 

4. The spectrum contains a soft excess with a thermal temperature of ~ 0.25 

keV, or a power law with an a of 4.5 ± 1.0, and a luminosity about 0.01 to 0.1 

the hard component, or^ broad iron i.-emission at ~ 1 keV with an EW of ~ 400 

eV. 

£)_ 4ll2129^^47 

The low state data come from a HEAO-1 A2 pointed observation in Dec 1978, 

and the high state data are from SSS/MPC and were taken six months later. 

2 

The MED gives a power law with a = 1.1 ± 0.2 for ^ ~ 0*8. 3 thermal 

2 

of temperature 7.3 ± 1.5 keV for " 1*3. For both spectral models a 425 ± 

250 eV EW iron line with a FWHM < 3.5 keV is required. The column density is 

< 5 X 10^^ H cm"^. The HED gives similar results, with a slightly 

better for the power law. The SSS spectrum, taken when the source was 

bright, indicates a flatter spectrum with a = 0.3 ± 0.3 or kT > 8.0 keV (x^ = 

2 

0.8 in both cases). The MPC, however, gives a = 1.0 ± 0.1 for x^. = 3.2, or kT 
= 8.9 1 1.0 keV for x^ = 0*7, in the 1.5 to 20 keV band. The inverted spectra 
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tor all three detectors can be compared in Figure 14. All the data are 
marginally consistent with a single thermal spectrum of ~ 8 keV. The high 
state flux in the 0.5 to 60 keV band is 3 x 10"^® ergs cm"^s“^, which gives a 
luminosity of 10^® ergs s"^ for a distance of 5 kpc. 

VIII. SPECTRAL INTERPRETATION 

All models for Cyg X-3, 4U1822-37 and 4U2129+47 must invoke the 
Comptoni-zation of photons from a central X-ray source to explain the X-ray 
light curves. The nature of this process in an astrophysical setting has been 
discussed by many authors (e.g. Illarionov and Sunyaev 1972; Felten and Rees 
1972; Shapiro, Lightman and Eardley 19/6; Katz 1975). 

£)_ 4U1822-37 and Cyg X-3 

The low state Cyg X-3 spectrum and the spectrum of 4U1822-37 are 
remarkably similar. For a scattering medium with a temperature which is well 
below the characteristic source temperature, the electron recoil effect will 
dominate the Compton scattering process, and a break will be seen in the 
spectrum at Ej^ = M^c'^/t (Illarionov and Sunyaev 19/2). The high energy 
cutoff in 4U1822-37 and the low state of Cyg X-3 is a direct measure of the 
optical depth to the central X-ray object, and gives, in both cases, t ~ 5. 

As T increases, the energy of the high energy cutoff will decrease, 
approaching the temperature of the corona. For i > 10 photons will become 
sufficiently themialized to form a Wien peak with the characteristic 
temperature of the corona. The high state of Cyg X-3 can be reconciled with 
an increase in t to > 10 that may be associated with an increase in the 
central source luminosity. The fact that there is a high energy tail 
Indicates a non-uniform or leaky scattering medium (which would be expected 
from an ADC, cf Figure 4). The blackbody spectrum gives a temperature 
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Of ~ 1.5 koY for the scattering medium. 

The 'act that both 4U1H22-3/ and Cyg X-3 display intense broad iron line 
emission at 6-/ keV further emphasises that these sources may be fundamentally 
similar. The EWs of order 1000 eV are at least a factor 2 greater than that 
seen from other binary X-ray sources. An optically thick ADC provides an 
ideal environment for the production of line emission, either from 
recombination in the corona (Ross, Weaver and McCray 1978; Ross 1979; Fabian 
and Ross 1981), or from the fluorescence of iron at the surface of the 
accretion disk (paper 1). Either process must account for the factor 2 
decrease in the EW during the high state of Cyg X-3. In the high state the 
reduced fluorescence efficiency, caused by the softer continuum above 7.1 keV, 
Is consistent with lower EW. However, the observed EWs requires, for solar 
abundances, that the flux at the surface of the disk is a factor 10 larger 
than that seen from our line of sight (paper 1); this is unlikely for the 
geometry discussed in V. Recombination in the corona may be preferable to 
fluorescence but the predicted EWs and spectra are critically dependent on the 
assumed density structure of the ADC, a highly uncertain parameter fer the 
Compton- thick case. 

The broadening of the line profiles is a direct consequence of the 
presence of the hot Compton- thick ADC. The iron line profile, including the 
energy centroid, conveys information conr»»rn1ng both the temperature and the 
optical depth of the scattering medium (Sunyaev and Titarchuk, 1980, ST; 
Pozdnyakov, Sobol, and Sunyaev 19/9, PSS). If the temperature of the corona 
is hotter or cooler than the line energy, the scattering will broaden the line 
profile and shift its centroid. The line energy and width of Cyg X-3 is the 
same in both the low and hign states, which infers kT^, ~ 2 keV for i between 5 
and 10 (cf Fir,ure 4 of ST and Figure 4 of PSS). This is in reasonable 
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agreement with the Wien peak temperature. The width of the line In 4U1822-37 
Is about a factor two larger and Infers kT^. ~ 4 keV. All these values are 
consistent with the coronal temperatures Inferred from the measured coronal 
**ad1 1 . 

In Paper 1 we suggested that the soft excesses from Cyg X-3 and 4U1822-37 
might be free-free emission from their coronae, with a kT^. of ~ 0.25 keV. 

Inis no longer appears to be viable, as the combination of the required t and 
such a low coronal temperature would cause the Iron line to be measurably 
shifted to lower energies. A conservative estimate of the emission measure EM 
of the corona can be obtained by assuming the vertical density structure is 
constant at n^^^^ and cuts off at 0.1 Rq. This gives a luminosity for the 
free-free emission with a temperature of a few keV about a factor 100 below 
that of the observed emission. 

The disk acts as a partition across the center of the corona which, 
because It absorbs at least ~ 301 of the scatterred radiation which Is then 
re-emitted as a far UV albedo (kT ~ 0.01 kT^), substantially amplifies (by two 
orders of magnitude) the number of photons In the corona. The Doppler effect 
will dominate the Compton scattering of the disk albedo, which cools the 
corona (e.g. Katz 1975; Shapiro, Lightman and Eardley 1976). The Comptonized 
spectrum Is relatively independent of the input spectrum and for low x It will 
be a power law fro«n the disk surface temperature up to ~ kT^,, above which It 
decays as a Wien tall. The energy Index a Is given by ST and Shapiro, 

Lightman and Eardley to be 
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An upper limit to t can be obtained from the point at which catastrophic 
cooling sets which from ST is v ~ 5. This predicts t = 8-11 for kT^. = 
2-5 keV, in reasonable agreement with the observed values. The soft excess 
from 4U1822-37 may be the high energy tail of the albedo component. The 
measured a for this component of 4.5+1 .0 is consistent with the predicted 
value of 5. Suci. a soft component from Cyg X-3 would be undetected because of 
the high column density to that source. 

The X-ray spectra of 4U1822-37 and Cyg X-3 are therefore consistent with 
being the result of the Comptonization of a hard X-ray spectrum in a hot 10^ K 
optically th'‘ck plasma, such as would be found in an ADC. Previously proposed 
models for Cyg X-3 cannot reproduce all the experimental details. Scattering 
in an optically thick stellar wind from the secondary or a shell surrounding 
the entire system has been proposed to account for the 4.8 hr modulation from 
Cyg X-3 (Davidsen and Ostriker 1974; Pringle 1974; Milgrom 1976). Eisner et 
al. (1980) and Ghosh et al. (1981) have considered these models in great 
detail and can reproduce the asymmetric light curve of Cyg X-3 with an 
eccentric orbit. Our data suggest the high and low state of Cyg X-3 involve 
changes in t which, for the stellar wind or shell models, implies that the 
high state light curve should be much less modulated, in contrast to the 
constant fractional modulation observed (Figure 9). Further, if the 
scattering is in a stellar wind then an increase in absorption at minimum 
light, an order of magnitude larger than our upper limit of < 1.5 x 10^^ H 
cm"^, is tv. be expected. Lastly, apsidal motion will cause the light curve to 
change and the lack of any such change over the past 5 years constrains the 
secondary to be a helium star. For an ad hoc ;he11 surrounding the system the 
helium star must be "dressed" to reproduce the duty cycle. There is too much 
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extiP tion to Cyg X-3 to check the nature of the secondary directly. For 
4U1822-37, however, a helium star cannot be reconciled with the observations. 
^ 4U2129^•47 

The relatively soft ~ 8 keV thermal character of the spectrum of 
41)2129+47 Is very reminiscent of that seen from the Sco X-1 /galactic bulge 
sources (e.g. Mason et al. 1976). In these sources deviations, from a pure 
bremsstrahlung model have been accounted for by Comptonizatlon In an optically 
thick X-ray emitting plasma (e.g. Lamb and Sanford 1979). The ADC model 
provides a natural environment for such a process to occur, and 4U21 29+47 may 
provide the key to understanding these enigmatic sources. The spectrum of 
4U2129+47 probably differs from that of 4U1822-37 and Cyg X-3 because of a 
generic difference In the nature of the central X-ray emitting region. For 
example the magnetic field on the accreting neutron star may be weak and 
allows the accretion disk to penetrate close to the neutron star surface. 
Because of this the optical depth of the corona In the central region may be 
very high. In which case the radiative transfer is dominated by bremsstrahlung 
losses from the corona Itself (Fabian, priv. comm.). 

IX LUMINOSITIES 

A possible distinguishing property between turbulent and evaporated 
coronae Is the source luminosity, I.e. an evaporated Compton thick corona 
should not be found in a source with L^^ < 10^^ ergs s~^. The observed 

scattered luminosity LqJjj will be lower than the true source luminosity L by 
the fractional solid angle subtended by the corona, sin (where Is 
the angle above the disk plane at which t = 1). From Figure 3 0 > 45°, so 

LqIjj ~ 0.7 L. This gives an average high state luminosity for Cyg X-3 
comparable to the Eddington limit for 1.4 M^ star, which suggests that this 
source contains an evaporated corona. 
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4U1822-37 Is at 4 = 357° and b * -11°, and only a lower limit to its 
i distance of 500 pc exists (Mason et al. 1980). If its luminosity is similar 

to the low state Cyg X-3 value then its distance is > 10 kpc and the source is 
a halo object somewhere on the other side of the galaxy (Paper 1). Recently 
Cowley. Crampton and Hutchings (1982) have reported that 4U1822-37 has 
peculiarly large space motion of -90 km s~^. This supports the idea that this 

I 

source is in the galactic halo. 

The distance to 4U21 29+47 is estimated to be between 2-5 kpc (McClintock, 
Remillard and Margon 1981), so that L ~ 10^® ergs s"^, which is too low for a 
Compton thick evaporated corona. Rather than invoking mechanical heating of 
the corona for 4U2129+47, it is possible that because the inclination of this 
system is > 82°, the central area of the disk and a large part of the corona 
are hidden behind the disk rim (analogous to a permanently turned off Her 
X-1). This may also explain why the optical light curve of 4U2129+47 is 
dominated by the X-ray-heated secondary, and not, as is the case for 
41)1822-37, the light from the accretion disk. Therefore all three sources may 
well contain evaporated ADC. 

X. SUMMARY AND CONCLUSION 
The principle results can be summarized as fellows: 

1) The X-ray emission from 4U2129+47 is diffused by an ADC and the X-ray 
source appears to be extended with a radius of 0.6 Rq, compared to ~ 0.3 Rq 
; for 4U1822-37. 

f 

i 

1 2) ADC are probably generated by the evaporation of hot material from the 

< 

surface of the accretion disk by the central X-ray source. An evaporated ADC 
will only become optically thick when approaches the Fddingt'm limit of a 
1.4 Mq neutron star. 

3) Secondary modulations of the light curves of 4U1822-37 and 4U2129+47 can 
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be understood as the partial occultatlon of the ADC by structure In the outer 
regions of the accretion disk arising from turbulence Induced by the Inflowing 
gas stream. 

4) Detailed modelling of the light curve of 4U1822-37 Indicates two major 
structures In the outer disk thickness, one at the Initial point of Impact of 
the gas stream with the disk and a second ~ 230° around from that point. 

5) This model can be extended to account for the 4.8 hr modulation of Cyg 
X-3. In this case the Inclination is such that the eclipse by the companion 
Is not seen. 

6) The X-ray spectra of Cyg X-3 and 4U1822-37 are very similar and can be 
qualitatively understood as resulting from the Comptonizatlon of a hard X-ray 
spectrum by an optically thick ADC. 

7) We set an upper limit to any phase related changes In absorption from Cyg 
X-3 of < 1.5 X 10^2 H cm'2. 

8) The spectrum of 4U2129+47 Is similar to that from the Sco X-l/galactic 
bulge sou»'ces and indicates that the central X-ray source In this system Is 
different from that In 4U1822-37 and Cyg X-3. 

9) The Inclinations of Cyg X-3, 4U1822-37 and 4U2129+47 are < 70°, ~ 75° 
and > 82° respectively. The X-ray heated central region of the accretion disk 
of 4U2129+47 Is probably hidden behind a thick rim around the periphery of the 
disk. 

The gamma ray emission from Cyg X-3 has often been cited as evidence for 
a young pulsar In this system (e.g. Mllgrom and Pines 1978). The ADC model 
does not specifically require an accreting neutron star behind the corona, 
however models for the gamma ray emission that Invoke Inverse Compton 
scattering of the X-ray emission on relativistic electrons can also be 
acconwodated In an ADC format (Meszaros, Meyer and Pringle 1977; Fabian, 
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Blandford and Hatchett 1977). Intense Infrared emission from Cyq X-3 is also, 
on occasions, modulated with the 4.8-hr X-ray period with a similar light 
curve, but of lower amplitude (Mason et al. 1976a). This suggests that the 
modulated IR emission comes from the ADC. The nature of this IR flux and the 
dramatic radio outbursts remains puzzling. 

Evaporated optically thick ADCs will only be found among the more 
luminous X-ray sources that are close to being Eddington-limited, which 
suggests that many of the bright galactic center sources may contain ADCs. 

Our conclusions provide independent justification for the Comptonization 
models which have been applied to these sources. Milgrom (1978) has suggested 
that the lack of eclipses in the galactic center sources is caused by a thick 
accretion disk shadowing the companion in these systems. Our modelling of 
accretion disks thickened by turbulence supports this idea, although the 
presence of ADCs will, in a few systems, cause partial eclipses either by the 
companion or the outer disk structure to be seen. These may, however, be 
among the fainter galactic center sources, because a large part of the corona 
could be hidden by the edge of the disk (similar to the situation that may 
exist in 4U2129+47). The recent detection of periodic absorption events from 
4U1915-05 (White and Swank 1902, Walter et al. 1982), which can be understood 
as occultations of a central point source by structure in outer regions of the 
disk similar to that founJ here, supports the above scenario. In this case 
the steady luminosity of the source is ~ 10^^ ergs s"^, which is below the 
threshold for the formation of an optically thick ADC. 4U1822-37, Cyg X-3 and 
4U2129+47 may represent the "missing link" in our understanding of the bright 
galactic center sources and their relation to the better-understood pulsating 
binary X-ray sources. Future investigations into both Compton scattering 
processes and the dynamics of disk accretion should concentrate on reproducing 
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the properties of these three systems. 
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TABLE OPTIMUM BULGE PROFILE FITS 


(a) 4U1822-37 

(b) Cyg X-3 

5 » 40° ± 2 

40° (F) 

Op - 190° i 15 

100° t 10 

Op « 235^ ± 20° 

290° ± 15 

0p2 * 2iO° ± 20 

160° t 15 

H “ 0.55 1 0.05 

P 

0.55(F) 

All uncertainties are from eye estimates. 

(F) Indicates the parameter was 

fixed. Hp Is the height of the secondary 

bulge relative to the primary. 
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FIGURE CAPTIONS 


Figure 1 - The HEAO-1 A2 scanning data (with background subtracted) from 

4U2129+47 folded at the optical period, with phase zero corresponding to 
optical minimum. The data Is repeated for half a cycle. 

Figure 2 - The SSS/MPC observation of 4U2129+47 folded at the optical 

period. The background has been removed and the plot repeated for half a 
cycle. 

Figure 3 - The optical depth through the corona as a function of viewing 
Inclination for various source luminosities, under the idealized 
condition that the radiation pressure Is zero. When t > 1 the source 
luminosity must be scaled down (by a factor we have not attempted to 
calculate here) In order to account for the enhanced photon density In 
the corona (Ross 1979). 

Figure 4 - Two idealized views of the possible appearance of the 4U18'Z2-37 

system. The cross sectional view of the disk assumes we are viewing 
at ~ 0.15. The hatched and dotted areas represent the optically thick 
and thin regions of the ADC, respectively. The binary phases are 
indicated for the top view. 

Figure 5 - Geometrical restrictions from the X-ray eclipses of 4U1822-37 
and 4U2129+47 with T7 set to 5 and 2, respectively. The area of 
uncertainty for 4U1022-37 contains both the statistical errors and the 
uncertainty in the mass of the neutron star. For 4U2129+47 only the 
statistical errors are included, and the uncertainty in the mass of the 
compact object will move the uncertainty region to the left or right by 
the Indicated amount. An additional uncertainty has been considered for 
4U1822-37 that pertains to the discussion of the possibility that there 
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is a thickening of the accretion disk at the edge caused by the inflowing 
material. The dashed lines give the region of uncertainty assuming that 
the disk thickness is much less than the coronal height. An accretion 
disk with appropriate "bulges" to generate the azimuthal profile for 
4Ulfl22-37 discussed in section V is indicated by the solid lines. 

Figure 6 - Four predicted light curve for the eclipse of an ADC by a bulge 
at the edge of the disk that illustrate the effect of changing the radius 
of the ADC Rj. for a constant disk radius R^j. The assumed bulge profile 
is shown at the bottom of the figure. The secondary eclipse has been 
eKcluded and the 4U1822-37 data shown for reference. 

Figure 7 - The same as Figure 6 with the addition of a second bulge. 

Figure 8 - Top two panels; the best fit to the HEAO-1 A2 41)1822-37 light 
curve for the indicated bulge profiles. Third panel down: the saine fit 

given in the second panel down compared to an Einstein observation 
one year later. Rottom panel; the predicted optical light curve 
compared to the observations from Mason et al . The dashed line 
represents the predicted eclipse assuming a uniformly illuminated disk 
interior and bulge. The solid line includes 40% of the light from a 
companion modulated by 100%, with minimum centered on the X-ray eclipse 
to approximate the effects of an X-ray heated companion. 

Figure 9 - The best fit eclipse profile for Cyg X-3 in its low state 

compared to Einstein MPC data (top panel). In the bottom panel the fit 
from the low state is compared to scanning data obtained in the high 
state. In both cases it is assumed that the inclination is too low for 
the eclipse by the companion to be seen. 

Figure 10 - Right panel: the allo»>.'ed region for the companion radius for 

Cyg X-3 as a function of inclination for M^^ = I.- deduced from the 
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lack of an eclipse. Left panel: the peak height of the bulge as a 

function of Inclination for both Cyg X-3 and 4U1822-37. 

Figure 11 - The deconvolved spectrum of Cyg X-3 for Its low and high states. 

Figure 12 - The PHA data from the minimum and maximum of the 4.8 hr cycle 

presented In the form of a ratio of the former to the latter. 

Figure 13 - Top Panel: the SSS PHA ratio defined In a way similar to that 
of Figure 12. Middle panel: the deconvolved spectrum for the 

two-component-plus-sul fur-1 ine fit. The dashed line Indicates the fit 
without the contribution from the s it component. Bottom panel: the 

best fit for a featureless power-1 aw-plus-absorptlon model. 

Figure 14 - The deconvolved spectra of 4U2129+47 for both the high and low 

states. 
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